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One  of  the  most  important  factors  affecting  the  attenuation  of  solar  radiation  reaching  the  earth’s  surface 
under  cloudless  sky  conditions  is  the  presence  of  aerosol  particles.  A  direct  measurement  of  the  aerosol 
transmittance  is  not  possible  due  to  the  strong  influence  of  the  other  atmospheric  components.  Thus 
the  extinction  caused  by  the  atmospheric  aerosol  can  be  calculated  only  indirectly  using  fundamental 
physical  techniques.  This  work  provides  an  overview  of  the  effect  of  aerosols  on  solar  radiation  budget 
by  considering  two  common  turbidity  parameters  including  the  Linke  turbidity  factor  TL  and  Angstrom 
turbidity  coefficients  y3  and  a.  Total  extinction  of  solar  radiation  due  to  the  absorption  and  the  scattering 
caused  by  the  atmospheric  aerosol  accounts  for  10-20%  for  zero  zenith  angle.  The  influence  of  aerosol  on 
radiation  passing  through  the  atmosphere  cannot  be  neglected,  especially  in  urban  or  industrialized  areas. 
The  attenuation  of  solar  radiation  through  a  real  atmosphere  versus  that  through  a  clean  dry  atmosphere 
gives  an  indication  of  the  atmospheric  turbidity. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

One  of  the  most  important  factors  affecting  the  amount  of  solar 
radiation  reaching  the  earth’s  surface  under  cloudless  sky  condi¬ 
tions  is  the  presence  of  aerosol  particles  (small  particles  either  in 
solid  state  or  in  liquid  state)  in  the  atmosphere  [1,2].  The  aerosol 
number  or  mass  concentration  is  key  parameters  in  the  studies  of 
atmospheric  radiative  effects  [3].  Aerosols  are  present  in  both  tro¬ 
posphere  and  stratosphere  and  mostly  throughout  the  atmospheric 
boundary  layer  at  number  concentrations  depending  upon  factors 
such  as  location,  atmospheric  conditions,  annual  and  diurnal  cycles 
and  presence  of  local  sources  [4].  Natural  aerosol  particles  range  in 
radius  from  1  to  105  nm;  very  small  particles  (called  Aitken  par¬ 
ticles)  from  1  to  lOnm;  and  large  particles  from  10* 1 2  to  103nm. 
Particles  in  the  103-105  nm  range  are  called  giant  particles  [5,6]. 
According  to  the  formation  of  the  aerosol  dispersions,  they  fall 
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into  two  categories,  viz.,  dispersed  aerosols  and  condensed  aerosol 
[2].  The  presence  of  aerosols  make  the  atmosphere  turbid  and  the 
determination  of  turbidity,  the  property  of  the  atmosphere  in  the 
presence  of  aerosols  is  of  importance  in  climatology,  in  pollution 
studies  and  for  solar  energy  utilization  [7].  Aerosols  influence  the 
solar  radiation  both  directly  and  indirectly  through  their  various 
sizes  and  thus  their  different  optical  and  physical  properties.  When 
aerosols  are  sufficiently  large  in  size,  they  scatter  and  absorb  sun 
light,  and  when  these  particles  are  small,  they  act  as  cloud  conden¬ 
sation  nuclei  and  aid  in  the  formation  of  clouds.  The  attenuation 
caused  by  absorption  and  scattering  by  aerosols,  is  known  as  atmo¬ 
spheric  turbidity.  The  assessments  of  radiative  effects  of  aerosols, 
especially  of  the  indirect  effect  on  the  radiative  balance  of  the  earth 
are  highly  uncertain.  Uncertainty  arises  in  part,  because  of  the  dif¬ 
ferences  in  the  assumed  aerosol  concentration  and  evaluation  of 
the  cloud  droplet  number  concentration  (CDNC)  [3].  Although  the 
amount  of  aerosols  presented  in  the  atmosphere  in  the  vertical 
direction  can  be  quantified  in  terms  of  the  number  of  particles 
per  unit  volume  [8]  or  their  mass  in  microgram  per  cubic  meter 
[7]  existing  techniques  for  determining  this  number  is  practically 
applicable  only  for  aerosols  presented  near  the  earth’s  surface  [8]. 
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Nomenclature 

P  the  Angstrom’s  turbidity  coefficient 

a  the  Angstrom  wavelength  exponent 

N  number  density  of  aerosols 

X  wavelength 

rx  aerosol  optical  depth  at  the  wavelength  X 

r  radius  of  aerosols 

v  fitting  parameter  characterizing  the  size  distribu¬ 

tion 

h  solar  elevation  angle 

Po  true  value  of  Angstrom’s  turbidity  coefficient 

a0  true  value  of  Angstrom  wavelength  exponent 

OG1  old  number  of  cutoff  filter  OG530,  opaque  up  to 
530  nm;  transparent  530-2800  nm 
RG2  old  number  of  cutoff  filter  RG630,  opaque  up  to 
630  nm;  transparent  630-2800  nm. 

RG8  old  number  of  cutoff  filter  RG695,  opaque  up  to 
695  nm;  transparent  695-2800  nm 
w  water  vapor  content 

Tl  the  Linke  turbidity  factor  for  an  air  mass  equal  to  2 

TLm  long-term  monthly  average  of  the  Linke  turbidity 
factor  for  an  air  mass  equal  to  2 
<5R(m)  the  integral  Rayleigh  optical  thickness 
<5cda  clear  and  dry  atmosphere  optical  thickness 
m  air  mass 

rL1P  Linke  turbidity  factor  corrected  by  Ineichen  and 

Perez 


However,  it  is  more  common  to  represent  the  amount  of  aerosols 
by  an  index  of  turbidity  [7]. 

Due  to  the  relationship  existing  between  aerosols  and  atten¬ 
uation  of  solar  radiation  reaching  the  earth’s  surface  [9]  several 
so-called  “turbidity  indices”  have  been  defined  during  the  past 
decades,  differing  slightly  in  physical  meaning.  The  basic  relation¬ 
ship  underlying  the  derivation  of  the  various  turbidity  indices  is 
the  Lambert-Bouguer-Beer  law  [10].  There  are  at  least  11  meth¬ 
ods  which  discussed  atmospheric  turbidity  in  the  literature  such 
as  Linke  and  Boda  [11],  Kastrov  [12],  Angstrom  [13],  Feussner  and 
Dubois  [14],  Schneider  [15],  Schuepp  [16],  Volz  [17],  Valko  [18], 
Unsworth-Monteith  [19],  Davies  and  Hay  [20],  illuminance  tur¬ 
bidity  [21].  Among  these  atmospheric  turbidity  coefficients,  the 
most  frequently  used  are  the  Angstrom  turbidity  coefficients  and 
the  Linke  turbidity  factor,  as  reported  by  Wen  and  Yeh  [22].  The 
Angstrom  exponent  is  very  popular  because  of  the  simplicity  of  the 
respective  equation,  because  it  enables  to  interpolate  or  to  extrap¬ 
olate  aerosol  optical  properties,  and  because  it  is  connected  to 
particle  microphysics  [23  ].  In  addition,  it  also  indicates  the  amounts 
of  aerosols  and  represents  the  combined  effects  of  both  scattering 
and  absorption  caused  by  aerosols.  It  also  indicates  the  amounts  of 
aerosols  [24].  The  Linke  turbidity  factor  TL  is  of  particular  impor¬ 
tance  for  the  calculation  of  beam  solar  and  diffuse  radiation  in  a 
determined  locality.  This  parameter  is  also  useful  for  the  predic¬ 
tion  of  the  availability  of  solar  radiation  and  day-light  illuminance 
under  cloudless  skies  [25]  and  it  is  a  key  input  to  several  mod¬ 
els  assessing  the  downwelling  irradiance  under  clear  skies  that  are 
used  by  several  communities  in  the  fields  of  renewable  energies, 
climatology,  agro-meteorology,  and  atmospheric  pollution  [26]. 

2.  Angstrom’s  turbidity  coefficients  P  and  a 

The  aerosol  optical  thickness  depends  not  only  on  aerosol  char¬ 
acteristics  (size  distribution,  refractive  index,  etc.)  but  also  on 
aerosol  total  loading.  Angstrom  suggested  a  single  formula  for 


aerosol  scattering  optical  thickness  evaluation  generally  known  as 
Angstrom’s  turbidity  formula  given  by  the  following  [27] 

rx  =  fa~a 

In  this  equation  p  is  called  the  turbidity  coefficient  and  a  the  wave¬ 
length  exponent. 

The  Angstrom  turbidity  coefficient  p  represents  the  amount  of 
aerosols  in  the  atmosphere  in  the  vertical  direction  [28 ].  The  param¬ 
eter  p ,  which  may  vary  from  0.0  to  0.5,  is  an  index  of  the  amount  of 
aerosols  present  in  a  vertical  column  of  the  atmosphere.  The  param¬ 
eter  a  is  a  reliable  index  of  the  size  distribution  of  these  aerosols  (at 
least  for  particles  of  radii  ranging  from  0.1  to  1.0  mm),  as  indicated 
by  the  nephelometer  measurements  of  Charlson  (1972)  cited  in  [1  ]. 
In  the  solar  spectrum,  the  Angstrom  exponent  a  is  a  good  indicator 
of  the  dominant  size  of  the  atmospheric  particles  [29].  Fluctuations 
of  a  reflect  aerosol  size  distribution  variations.  The  maximum  value 
of  a  (equal  to  4)  corresponds  to  molecular  extinction.  Values  near 
zero  (or  even  negative)  correspond  to  coarse-mode  aerosols  (sea 
spray  and  desert  dust)  indicating  a  non-aerosol  optical  depth  (r) 
wavelength  dependence,  while  values  of  a  above  1.5  indicate  sig¬ 
nificant  presence  of  fine-mode  particles  (mainly  smoke  or  urban 
aerosols).  The  average  a  over  the  whole  spectrum  between  290  and 
4000  nm  has  a  value  between  0.9  and  2.0  with  the  most  frequent 
observations  around  a  =  1.3  [30]. 

Several  methods  may  be  used  to  determine  the  coefficient  P 
and  even  the  coefficient  a.  Some  of  them  only  require  knowledge 
of  the  direct  solar  total  irradiance  data,  while  the  other  need  one  or 
two  spectral  irradiance  measurements  [1  ].  They  are  formulated  by 
Louche’s  et  al.  [31],  Gueymard  and  Vignola  [32],  Pinazo  et  al.  [33], 
Grenier  et  al.  [34],  Method  U  as  stated  by  Casiniere  et  al.  [1],  Volz 
method  [17],  direct  method  and  linear  fitting,  as  cited  by  Cachorro 
et  al.  [35],  spectral  window  method,  as  stated  by  Kaskaoutis  et  al. 
[36].  Junge  [37]  was  the  first  to  explore  a  relationship  between  the 
Angstrom  exponent  a  and  the  aerosol  size  distribution  by  assuming 
a  power-law  relationship  for  the  number  density  (N)  of  aerosols  as 
function  of  radius  (r).  Additionally,  using  Mie  theory  for  spherical 
particles  in  the  range  1 00  nm  <  r  <  1 000  nm,  Junge  was  able  to  show 
that  a  «  v  -  2  for  nonabsorbing  aerosols  with  a  >  1  as  reported  by 
Kaskaoutis  et  al.  [38].  According  to  EI-Hussainy  and  Omran  [39] 
Angstrom  [40]  discussed  a  rather  simple  method  to  determine  the 
true  values  of  both  Angstrom  parameters  p0  and  uq  based  on  the 
assumption  that  the  parameter  p  is  constant  for  all  points  of  the 
spectrum  provided  it  is  computed  by  using  the  real  value  of  a. 
Various  studies  [41-43]  as  reported  by  [37]  have  shown  that  the 
Angstrom  formula  is  a  special  case  of  a  more  complicated  law  valid 
for  a  limited  range  of  values  on  the  diameter  of  particles  and  for 
a  limited  interval  of  wavelengths.  A  method  of  determining  the 
true  turbidity  coefficient  p0,  corresponding  to  the  real  value  of  a , 
was  proposed  by  Herovanu  [44]  and  simplified  by  Angstrom  [45] 
as  reported  by  Joseph  and  Manes  [10].  This  method  is  based  on  the 
assumption  that  P  and  a  of  Angstrom’s  empirical  aerosol  extinc¬ 
tion  law  can  be  determined  by  considering  the  intensities  measured 
within  two  spectral  ranges,  0-525  nm  and  525-630  nm,  with  OG1 
and  RG2  cutoff  filters.  Additionally,  Rodhe  [46]  has  proposed  a  com¬ 
putational  procedure  for  Angstrom  turbidity  coefficient  p.  Volz  [47] 
proposed  it  is  possible  to  determine  the  values  of  p  and  a  simulta¬ 
neously  from  aerosol  transmittances  at  two  wavelengths  using  sun 
photometers  as  cited  in  Utrillas  [48].  According  to  Casiniere  et  al. 
[1  ],  Dogniaux  [49]  developed  the  relationship,  as  reported  by  Page 
[50]  from  a  statistical  analysis  of  the  many  p  values  determined 
in  meteorological  stations  (by  means  of  pyrheliometers  fitted  with 
standard  color  filters  such  as  OG1,  RG2,  and  RG8). 

King  and  Byrne  [51]  showed  that  the  size  distribution  of  aerosols 
does  not  typically  follow  the  Junge  law  nor  does  it  take  values  from 
zero  to  unity  as  reported  by  Kaskaoutis  et  al.  [52].  According  to 
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Kaskaoutis  et  al.  [37]  King  and  Byrne  [51  ]  have  shown  that  aerosol 
size  distributions  do  not  have  radii  extending  from  zero  to  infinity. 
Louche  et  al.  [31  ]  reported  King  and  Buckius,  based  on  the  work  of 
others,  have  given  a  simple  expression  relating  visibility  to  a  and 
p.  With  the  help  of  the  Bird  and  Hulstrom  [53,54]  parameteriza¬ 
tion  of  the  direct  solar  total  irradiance,  Louche  et  al.  [31  ]  described 
a  simple  approach  to  determine  Angstrom  turbidity  coefficient  (3 
under  cloudless  skies.  Melice  et  al.  [55]  calculated  turbidity  coef¬ 
ficients  in  the  Sahara  of  Tozeur  in  southern  Tunisia.  They  found 
high  turbidity  values  in  spring  and  summer.  Since  these  turbidity 
indices  differed  slightly  in  their  meaning,  correlations  and  relation¬ 
ships  were  established  between  them  as  stated  by  EI-Hussainy  and 
Omran  [39]. 

Casiniere  et  al.  [  1  ]  reported  Nicholls  [56]  simultaneously  derived 
the  dependence  of  the  parameter  a  on  the  wavelength.  Bird  [57] 
observed  that  turbidity  versus  wavelength  often  exhibits  a  curve¬ 
like  relationship  on  a  log-log  plot.  Different  procedures  can  be  used 
as  stated  by  Cachorro  [35]  and  Kaskaoutis  et  al.  [52]  for  the  purpose 
of  determining  Angstrom  parameters  namely,  Volz  method,  direct 
method,  linear  fitting  and  spectral  window  method.  EI-Hussainy 
and  Omran  [39]  reported  Louche  et  al.  [31  ]  and  Fox  [58]  estimated 
Angstrom  turbidity  parameters  using  measurements  of  total  spec¬ 
trum  direct  irradiance.  O’Neill  and  Royer  [60]  analyzed  a  computed 
from  five  VIS  and  NIR  wavelengths  (380-862  nm),  and  found  a 
relationship  between  the  effective  radius  in  the  aerosol  size  dis¬ 
tribution  and  the  value  of  a  as  cited  by  Kaskaoutis  et  al.  [59].  As 
reported  by  Basart  et  al.  [29],  Kaufman  [61  ]  pointed  out  that  nega¬ 
tive  values  of  the  difference  8a  =  <2(440,61 3)  -  of(61 3,1 003)  indicate 
the  dominance  of  fine  mode  aerosols,  while  positive  differences 
reflect  the  effect  of  two  separate  particle  modes. 

Gueymard  [62,63]  expressed  the  turbidity  as  Bird  [57]  did,  two 
different  values  for  a:  a\  for  X  <  500  nm  and  a2  for  X  >  500  nm  as 
reported  by  Casiniere  et  al.  [1  ].  As  cited  in  Janjai  [24]  Kaufman  et  al. 
[64]  argue  the  need  to  measure  aerosol  characteristics  at  different 
locations  and  in  various  meteorological  conditions  to  fully  assess 
regional  aerosol  properties.  Grenier  et  al.  [34]  presented  three 
linear  equations  valid  for  three  ranges  of  columnar  water  vapor 
content  w  for  calculating  the  Angstrom  turbidity  coefficient  ft  from 
Tl  as  reported  by  Kasten  [65].  Pinazo  et  al.  [33]  propounded  new 
method  for  calculating  the  Angstrom  turbidity  coefficient  based  on 
the  ratio  of  direct  solar  radiation  to  global  solar  radiation  on  a  hor¬ 
izontal  surface  and  on  the  “C”  model  of  Iqbal.  Utrillas  [47]  reported 
Bokoye  et  al.  [66]  have  shown  that  the  optimum  wavelength  pair 
in  the  290-900  nm  interval  is  400  nm  and  750  nm  when  using  the 
Volz  method.  Gueymard  [32]  has  proposed  a  semi-physical  method 
to  evaluate  turbidity  from  broadband  irradiance  measurements 
and  other  atmospheric  parameters.  This  method  requires  cloud¬ 
less  skies  and  very  accurate  measurements  of  global  (or  diffuse) 
and  direct  radiation.  Eck  et  al.  [67]  showed  how,  in  the  wavelength 
range  of  340-870  nm,  a  can  increase  by  a  factor  of  3-5  as  wave¬ 
length  increases  for  biomass  burning  and  urban  aerosols,  while 
remaining  constant  or  decreasing  in  the  presence  of  mineral  dust. 

Malik  [68]  modified  Pinazo  et  al.’s  method  and  developed  two 
new  correlations  that  link  Pinazo’s  method  with  model  ‘C  by  Iqbal. 
The  modified  model  was  used  to  compute  Angstrom  turbidity 
coefficient  for  Brunei  Darussalam  using  measured  data  on  solar 
radiation  for  a  period  of  5  years  (1990-1994).  Adeyewa  and  Balo- 
gun  [69]  reported  Adeyefa  et  al.  [70]  also  examined  spectral  and 
broad-band  solar  radiation  measurements  made  non-continuously 
between  1990  and  1994  at  Abisko.  The  study  showed  that  there 
was  an  increase  in  atmospheric  turbidity  due  to  the  huge  amounts 
of  volcanic  aerosols  injected  into  the  stratosphere  by  the  volcanic 
eruptions  of  Mount  Pinatubo  in  June  1991.  Utrillas  [48]  devel¬ 
oped  a  new  method  aimed  at  determining  the  Angstrom  turbidity 
coefficients  from  measurements  of  broadband  direct  irradiance 
with  filters.  The  results  obtained  using  this  method  have  been 


compared  with  those  obtained  using  the  Angstrom  method,  and 
both  have  been  compared  with  the  values  obtained  from  1000- 
nm  wavelength  spectral  irradiance  measurements.  O’Neill  et  al. 
[71  ]  demonstrated  that  an  Angstrom  exponent-based  separation  of 
coarse  from  fine  mode  contribution  to  r  is  feasible  in  part  because  of 
the  coarse  mode  r  spectral  variation  being  approximately  neutral. 
Schuster  et  al.  [72]  addressed  the  link  between  Angstrom  expo¬ 
nent  curvature  and  the  ratio  between  fine  and  total  aerosol  volume. 
Recently,  Gobbi  et  al.  [73]  introduced  a  straight-forward  graphi¬ 
cal  frame-work  that  allows  to  discriminate  different  aerosol  types 
based  on  aerosol  spectral  measurements  by  sunphotometers  which 
can  be  characterized  by  three  independent  pieces  of  information: 
r,  a  and  the  spectral  curvature  of  a  (<$o0.  Plotting  data  in  this  space 
allows  for  inference  of  the  aerosol  fine  mode  size  and  fractional 
contribution  to  total  r. 

Tables  1-4  describe  summarizing  studies  which  have  been  con¬ 
ducted  at  different  regions  by  using  various  methods  that  are 
mentioned  above  to  determine  Angstrom  turbidity  coefficients. 

3.  The  Linke  Turbidity  coefficient  (Tl) 

The  Linke  turbidity  factor  TL  is  an  indicator  of  the  number  of 
clean  dry  atmospheres  that  would  be  necessary  to  produce  the 
attenuation  of  the  extraterrestrial  radiation  that  is  produced  by 
the  real  atmosphere.  Typical  values  range  from  1  to  10  [28].  TL 
can  be  obtained  directly  from  observations  performed  during  very 
clear  sky  periods,  but  this  kind  of  experimental  data  is  rarely  avail¬ 
able,  thus  Tl  is  generally  an  estimated  parameter.  Also,  time  series 
of  radiation  data  are  generally  too  short  to  allow  estimation  on  a 
daily  basis.  Fortunately,  long-term  monthly  average  values,  TLm, 
are  sufficient  for  most  applications  [26].  Feussner  and  Dubois  [14] 
published  a  series  of  spectral  data  tables  enabling  the  calculation 
of  SR(m )  where  both  molecular  scattering  and  absorption  by  the 
stratospheric  ozone  layer  are  taken  into  account.  Additionally,  Dog- 
niaux  [89]  derived  empirical  relation  from  extensive  experimental 
campaigns,  illustrating  the  observed  variation  of  Linke’s  turbidity 
factor  with  solar  elevation  angle,  h  (°),  atmospheric  water  vapor 
content  w  (cm)  and  Angstrom’s  turbidity  coefficient  (3  as  reported 
by  Molineaux  et  al.  [90].  According  to  Kasten  [65]  in  order  to  facili¬ 
tate  the  determination  of  T(m)  from  pyrheliometric  measurements, 
Kasten  [91]  presented  a  simple  analytic  formula  for  <$R(m),  based 
on  numerical  values  published  by  Feussner  and  Dubois  [14]  and 
internationally  recommended  by  CSAGI  [92].  Based  on  their  new 
values  of  5R(m)  Louche  et  al.  proposed  a  new  analytic  formula  for 
the  function  <5R(m)  and  demonstrated  that  the  new  <5R(m)  values 
are  definitely  higher  than  the  old  CSAGI  values  used  by  Kasten, 
particularly  at  small  relative  optical  air  masses  m. 

Ineichen  and  Perez  [93]  reported  that  Kasten  [91  ]  fitted  the  fol¬ 
lowing  equation  to  those  tables  published  by  Feussner  and  Dubois 
[14] 

Scda  =  (9-4  +  0.9  x  AM)~^ 

which  is  known  as  Kasten’s  pyrheliometric  formula.  Additionally, 
in  an  attempt  to  improve  the  formulation  proposed  by  Kasten  [91  ], 
Louche  et  al.  [94]  and  Grenier  et  al.  [34,95]  added  absorption  by  the 
permanent  gaseous  constituents  to  the  definition  of  <5R(m)  (these 
gases  are  considered  uniformly  mixed  and  invariable  in  both  a  clean 
dry  atmosphere  and  a  turbid  atmosphere).  As  cited  in  Navvab  et  al. 
[21],  Valko  [96]  has  suggested  empirical  fits  of  the  Linke  turbid¬ 
ity  factor  (Tl)  as  a  function  of  aerosol  and  water  vapor  content 
and  relative  air  mass.  Grenier  et  al.  [34]  reported  in  1986  that  a 
determination  of  <5R(m)  based  on  more  accurate  values  of  spec¬ 
tral  extraterrestrial  solar  irradiance  and  extinction  coefficients  of 
the  various  attenuators  was  carried  out  by  Louche  who  proposed 
the  following  algorithm  to  evaluate  the  optical  thickness  of  the  CD 
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Studies  of  atmospheric  turbidity  in  Europe  using  Angstrom  formula. 


4267 


Author 

Kaskaoutis  et  al.  [59] 


Kaskaoutis  and  Kambezidis  [36] 


Cucumo  et  al.  [74] 


Wagner  and  Silva  [75] 


Year  Method  Area  Results  and  conclusion 


2006  Based  on  ground  spectroradiometric 

measurements  taken  in  the 
atmosphere  of  Athens  during  May 
1995  Angstrom  exponent  values,  a, 
have  been  derived  at  five  narrow 
spectral  bands,  340-380,  380-440, 
440-500,  500-670  and  670-870  nm 
using  the  Volz  method  as  well  as  in  the 
whole  spectrum  (340-870  nm)  using 
both  Volz  and  least-squares  methods 


2008  Used  spectral  direct-beam  irradiance 

measurements  at  five  narrow  spectral 
intervals  (i.e.  340-380  nm, 

380-440  nm,  440-500  nm, 

500-675  nm,  675-870  nm)  as  well  as 
the  whole  spectrum,  320-1000  nm 
taken  in  Athens  during  May  1995,  to 
investigate  the  ability  of  the  different 
techniques  mainly;  the  Volz  method 
(VM),  the  direct  method  (DM)  and  the 
least-squares  method  (LM)  to  derive 
similar  Angstrom  turbidity  coefficients 
and  their  dependence  on  the  spectral 
range  used  for  their  determination 


1999  Carried  out  experimental 

determination  and  estimation  of  the 
Linke  turbidity  factor  and  Angstrom 
turbidity  coefficient  using 
experimental  data  of  beam  solar 
radiation  taken  at  Arcavacata  di  Rende, 
at  the  University  of  Calabria  from 
August  1995  to  August  1996,  and  at 
Casaccia  (Rome)  for  the  period  May 
1992-March  1995,  and  a  method 
introduced  recently  by  Pinazo  et  al.  for 
A  and  Dogniaux’s  formula  for  TL 


2008  Performed  a  simulation  study  in  order 

to  show  the  influence  of  the  aerosol 
optical  depth  (t)  distribution  together 
with  the  corresponding  error 
distribution  on  the  resulting  Angstrom 
exponent  (AE)  distribution 


Greece  1 .  a  values  depend  strongly  on  the 

wavelength  interval  used  for  their 
determination.  Very  high  a  values 
(2.83  ±  0.66)  were  derived  in  the 
340-380  nm  band,  while  as  the 
wavelength  interval  used  for  a 
determination  shifts  toward  longer 
wavelengths,  lower  a  values  were 
computed,  ( a  =  0.53  ±  0.39)  in  the 
670-870  nm 

2.  The  computed  a  values  in  the  whole 
spectrum  are  positively  correlated 
with  the  other  a  determined  in  the 
narrow  spectral  intervals 

3.  The  two  methods,  Volz  and  least 
squares  fit,  result  in  the  same  a  values 
in  the  340-870  nm  region 

Greece  1 .  The  various  techniques  lead  to 

different  Angstrom  turbidity 
coefficients  especially  if  narrow 
spectral  bands  at  the  shorter 
wavelengths  are  used 

2.  The  Angstrom  turbidity  coefficients 
derived  by  any  of  the  three  methods  at 
short  wavelengths  are  not 
representative  of  the  whole  spectrum 

3.  As  regards  the  whole  spectrum,  both 
LM  and  DM  estimated  quite  similar  /3 
and  a  values,  with  the  standard  errors 
of  each  fit  to  be  almost  the  same 

4.  On  the  other  hand,  when  narrow 
spectral  intervals  are  used  for  their 
determination,  the  errors  in  A  and  a 
increase 

Italy  1 .  The  Dogniaux  formula  allowed  the 

calculation  of  the  Linke  turbidity  factor 
with  average  errors  of  about  20%  and 
quadratic  errors  of  about  30%  assuming 
a  suitable  constant  value  of  the 
Angstrom  coefficient 
2.  The  Angstrom  turbidity  coefficient 
can  be  calculated  in  a  simple  way, 
using  the  Katz  formula,  with  an 
average  error  of  about  10%  and  an 
average  quadratic  error  of  30%  if  the 
experimental  value  of  the  Linke 
turbidity  factor  is  inserted  into  the 
Katz  formula 

Portugal  1 .  When  the  r  frequency  distribution  is 

wide  (e.g.  sigma  =  2)  then  the  resulting 
Angstrom  exponent  frequency 
distribution  is  narrower  in  comparison 
with  a  narrow  AOD-distribution  (e.g. 
sigma  =  1.1) 

2.  If  the  mean  z  is  high  (e.g.  0.4)  the 
Angstrom  exponent  frequency 
distribution  is  narrower  as  for  low 
turbidity  values  (e.g.  0.06)  in  case  of 
identical  absolute  errors 

3.  If  the  relative  errors  at  both 
wavelengths  are  equal,  or  in  a  practical 
sense  similar,  then  the  peak  or 
maximum  of  the  Angstrom  exponent 
distribution  reflects  the  true  value, 
otherwise  a  shift  either  to  smaller  or  to 
higher  Angstrom  exponent  values  will 
occur 

4.  Non  symmetrical  errors  result  in  a 
shift  of  the  maximum  of  the  Angstrom 
exponent  distribution 

5.  The  Angstrom  exponent  distribution 
is  less  symmetric  (higher  skewness) 
when  the  z  distribution  is  narrow  and 
more  symmetric  for  a  wide  z 
distribution 
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Table  1  ( Continued ) 

Author  Year  Method  Area 

Kaskaoutis  and  Kambezidis  [76]  2006  Compared  the  aerosol  optical  depths  (r)  Greece 

derived  by  the  Angstrom’s  formula  by 
means  of  direct-beam  spectra  measurement 
obtained  by  a  LICOR  spectroradiometer  in 
the  300-1 1 00  nm  wavelength  range  at 
several  sites  in  the  Athens  basin  under 
cloudless  skies  and  different  atmospheric 
conditions  during  May  1995  with  the 
spectral  aerosol  optical  depths  (SAODs) 
retrieved  through  the  SMARTS  model  and, 
therefore,  to  check  the  validity  of  the  former 
From  the  available  SAODs  derived  from 
SMARTS  model,  the  AODs  at  two  selected 
wavelengths  (500  and  1000  nm)  were  used 
for  further  analysis 


Polo  et  al.  [77]  2009  Determined  the  Angstrom  parameters  and  Spain 

total  ozone  column  obtained  by  means  of  a 
high-resolution  spectroradiometer  of 
selected  clear  sky  days  during  one  year  for  a 
characteristic  semi-desertic  site  in  the 
south-east  of  Spain 


Pedros  at  el.  [78]  1999  Determined  and  analyzed  the  Angstrom  Spain 

turbidity  coefficient,  the  Linke  turbidity 
factor,  and  the  Unsworth-Monteith 
coefficient  based  on  measurements  of 
normal  direct  irradiance  and  global 
horizontal  irradiance  taken  in  Valencia, 

Spain,  between  January  1990  and  December 
1996 


Canada  et  al.  [79]  1993  Employed  direct  radiation  measurements  Spain 

performed  in  Valencia  to  derive  the  monthly 
values  of  Angstrom’s  turbidity  coefficient,  A, 
under  cloudless  skies  by  means  of  Louche 
et  al.  method  and  then  compared  with 
models  reported  for  Avignon,  Ajaccio  and 
Dhahran 


Cachorro  et  al.  [80]  2001  Used  two  methods  namely,  the  direct  Spain 

method  and  the  window  method  and 
spectroradiometric  direct  irradiance 
measurements  at  four  non-absorbing 
spectral  windows:  A  (370-490  nm),  B 
(749-754  nm),  C  (776-782  nm),  and  D 
(862-872  nm)  under  clear  sky  conditions  at 
two  different  climate  stations  in  Spain,  in  a 
study  of  the  variation  in  the  Angstrom 
turbidity  parameter  a  and  its  dependence 
on  the  spectral  range  used  in  its 
determination 


Results  and  conclusion 

1.  The  Angstrom’s  formula  is  capable  of 
estimating  SAOD  values 

2.  The  derivation  of  the  r  values  obtained 
from  the  estimation  of  a  and  /3  values  was 
accurate  as  it  correlated  well  with  the 
SMARTS-estimated  SAOD  values 

3.  The  correlation  coefficient  between  SAOD 
and  z  was  always  above  0.7  and  in  the 
majority  of  the  cases  its  values  were  above 
0.9 

4.  The  correlation  coefficient  was  found  to 
depend  strongly  on  the  Angstrom’s  turbidity 
parameter  A  and  to  increase  as  /3  increases 

5.  At  500  nm  the  z  estimation  overestimated 
the  derived  SAODs  from  the  model  resulting 
in  negative  difference  values 

6.  At  1000  nm,  the  respective  differences 
were  higher  than  those  at  500  nm 

1.  The  turbidity  and  wavelength  exponent 
were  coherent  with  those  reported  by 
different  authors  for  desertic  area 

2.  A  seasonal  behavior  of  relative  low 
turbidity  values  around  0.02-0.1  was 
identified  during  the  autumn  and  winter 
days,  and  higher  values  in  the  range  of 
0.2-0. 6  were  found  for  spring  and  summer 

3.  Some  events  of  high  turbidity  (/I  around 
0.3-0.4)  with  low  wavelength  exponent 
(around  0.4-0. 6  p,m)  have  been  observed 
eventually  for  several  days  in  spring  and 
summer 

4.  Good  agreement  between  total  ozone 
column  values  determined  experimentally 
with  the  daily  evolution  reported  by  TOMS 

1.  The  values  of  all  the  considered  turbidity 
coefficients  have  a  minimum  in  the  winter 
and  grow  progressively,  reaching  a 
maximum  in  the  summer 

2.  The  summer  values  are  at  least  twice  as 
high  as  the  winter  ones 

3.  The  morning  turbidity  values  are  higher 
than  the  afternoon  ones,  and  all  three 
turbidity  coefficients  vary  in  a  similar 
manner  throughout  the  day 

4.  The  correlation  between  the  coefficients 
is  very  good  with  correlation  coefficients 
near  to  unity 

1.  The  turbidity  coefficient  of  Valencia 
showed  higher  and  lower  values  of  A  than 
Avignon  and  Ajaccio  in  summer  and  winter 
respectively.  This  is  mainly  due  to  the  effect 
of  the  air  mass  origin,  the  air  temperature 
and  the  water  vapor  content 

2.  The  variation  in  the  monthly  average 
values  at  Valencia  was  of  a  similar  trend  to 
that  of  Avignon,  Ajaccio  and  Dhahran 

1.  The  a-values  determined  in  the  UV  were 
very  different  than  those  determined  in  the 
VIS  and  VIS-NIR  ranges 

2.  The  low  values  determined  in  the  VIS-NIR 
interval  contrast  with  the  high  values  given 
by  the  350-400  nm  interval 

3.  The  differences  between  experimental 
and  estimated  values  were  as: 

•  For  the  340-400  nm  range  was  1% 

(s.d.  =  7%)  and  -2%  (s.d.  =  22%)  for  the 
350-400  nm  range,  respectively 

•  For  the  VIS  and  VIS-NIR  spectral  ranges 
gave  similar  results 

•  The  VIS  range  had  a  mean  value  of  26% 

(s.d.  =  15%)  and  the  VIS-NIR  is  32% 

(s.d.  =  16%) 

4.  These  results  seem  to  indicate  that  the  VIS 
and  VIS-NIR  ranges  are  not  useful 
predictors  of  the  AOD  in  the  UV. 
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Table  2 

Studies  of  atmospheric  turbidity  in  Asia  using  Angstrom  formula. 
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Author 


Year  Method 


Area  Results  and  conclusion 


Janjai  et  al.  [8] 


Hussain  et  al.  [7] 


Janjai  s.  et  al.  [24] 


Li  and  Lam  [28] 


Ganesh  et  al.  [81  ] 


Balakrishnalah  et  al.  [82] 


2003  Determined  Angstrom’s  turbidity  coefficient,  /l,  for 
53  meteorological  stations  covering  Thailand  by 
using  three  different  methods  mainly:  (a)  a 
calculation  of  /3  based  on  narrow-band  spectral 
irradiance;  (b)  a  derivation  of  /3  from  broad-band 
direct  irradiance;  and  (c)  an  estimation  of  /3  from 
visibilities 


2000  Determined  atmospheric  turbidity  parameters 
following  Angstrom’s  method  with  the  help  of  a 
pyrheliometer  and  cut-off  glass  filters  for  several 
months  of  1990, 1991, 1995  and  1997  for  Dhaka 
(latitude  23.7°N)  and  for  the  months  of  March  and 
April  1990  for  two  rural  locations,  Sripur  (24.1 1°N) 
and  Haripur  (26.03°N),  respectively 

2009  Employed  Bouguer’s  law  to  examine  aerosol 

optical  depth  and  Angstrom  coefficients  for  three 
sites  in  Bangkok  and  suburbs  under  cloudless 
conditions 

Three  types  of  Sun  photometers  have  been  used  to 
measure  direct  normal  spectral  irradiance  in  six 
wave-length  bands  at  a  frequency  of  1  min  at  these 
sites  for  a  period  of  2  years  (2004-2005);  Yankee 
MFRSR  for  Nakhon  Pathom  (NP)  at  wavelengths 
413,  500,  613,  671,  863  and  940  nm,  EKO  MS-110 
for  Phatum  Thani  (AIT)  at  wavelengths  368,  500, 
675,  778  and  863  nm  and  CIMEL  CE  317  for 
Bangkok  (BK)  at  wavelengths  380,  440,  500,  670, 
870  and  937  nm. 

2002  Used  three  methods  namely,  Louche’s  model  and 
Pinazo’s  model  and  visibility  method  to  determine 
Angstrom’s  turbidity  coefficient  (i.e.  /3 lou ,  /W  and 
/3Vis)  based  on  the  9-year  (1991-1999)  horizontal 
global  and  diffuse  solar  irradiances  measured  at 
the  City  University  of  Hong  Kong.  With  the  three 
sets  of  Angstrom’s  turbidity  data,  the  Linke 
turbidity  factors  (i.e.  TLo u,  TPin  and  Tv[s)  and  the 
original  Linke’s  expression  (TLin)  have  also  been 
computed  using  Dogniaux’s  formula.  In  addition, 
the  Linke  turbidity  also  determined  based  on  the 
original  Linke’s  expression  (TLin) 


201 1  Angstrom’s  turbidity  parameters,  a  and  |3  are 

calculated  and  analysis  on  daily,  monthly,  seasonal 
and  annual  using  the  solar  radiation  data  collected 
at  five  optical  channels  at  a  continental  station 
Mysore,  India  during  December  2003-June  2006 
performed  by  a  portable  sunphtometer 
MICROTOPS  II 


201 1  Used  Multi-Wavelength  solar  Radiometer  (MWR) 
measurements  preformed  during  2007-2008  over 
tropical  semiarid  site  of  Anantapur,  India  for 
determining  the  temporal  variability  of  aerosol 
characteristics 


Thailand  1 .  For  most  stations  in  the  north,  the  northeast  and 

the  central  region,  values  of  /3  were  relatively  high 
in  dry  season  (November-April)  and  low  in  the  wet 
season  (May-October) 

2.  For  the  dry  season,  it  was  inferred  that  the 
northeast  monsoon  together  with  the  high 
convection  in  the  summer  was  responsible  for  the 
higher  values  of  /3  for  these  regions 

3.  For  the  lower  values  of  /3  in  the  wet  season,  it  is 
due  to  the  effects  of  rains  which  remove  aerosols 
from  the  atmosphere 

4.  For  the  south,  values  of  /3  are  relatively  low,  and 
remain  nearly  constant  for  the  whole  year.  It  is 
presumably  caused  by  rains  from  the  two 
monsoons,  which  constantly  wash  through  the 
atmosphere  almost  all  year  round 

Bangladesh  1 .  Considerable  variation  of  Angstrom’s  turbidity 
parameters  for  Dhaka  over  the  year  was  observed 
with  a  maximum  value  in  March 

2.  The  value  for  Sripur  for  the  month  of  March  was 
somewhat  lower 

3.  The  value  for  Haripur  for  the  month  of  April  was 
not  much  different  from  that  of  Dhaka 

Thailand  1 .  Aerosol  optical  depth  showed  a  marked  change 

for  all  bands.  They  were  lowest  in 
October/November  and  gradually  increase  during 
the  dry  season,  reaching  maximum  values  in  late 
March  and  early  April 

2.  The  Angstrom  exponent  a  also  undergoes  a 
seasonal  variation,  with  maximum  values  in  the 
dry  season  and  minimum  values  at  the  height  of 
the  wet  season 

3.  There  was  much  more  scatter  in  the  seasonal 
pattern  of  the  turbidity  coefficient  /l,  with 
maximum  and  minimum  peaks  occurring  at 
various  times  of  the  year 

Hong  Kong  1 .  The  annual  averages  of  the  Angstrom  turbidity 
coefficients  based  on  the  Louche,  Pinazo  and 
visibility  methods  were  0.132,  0.165  and  0.156, 
respectively 

2.  In  general,  /3P in  produces  the  highest  monthly 
average  values  followed  by  /3Vis  and  then  /3L ou 

3.  The  ILin  has  the  lowest  monthly  average  values 
among  the  four  models,  ranging  from  3.7  to  5.26 

4.  For  the  other  three  models,  the  minimum 
monthly  average  occurs  at  TLou  in  December  at 
4.01  while  the  monthly  peak  appears  at  TPin  in  May 
at  6.72 

5.  It  is  envisaged  that  in  Hong  Kong  the  clear  sky 
conditions  can  be  defined  as  between  turbid  and 
clear 

India  1 .  The  observations  show  that  a  and  |3  vary 

throughout  on  individual  day  because  of  changes 
in  the  atmospheric  meteorological  parameters 

2.  It  is  observed  that  [3  is  highest  during  summer 
and  lowest  during  winter 

3.  An  anti-correlation  between  a  and  |3  is  observed 
throughout  the  day  during  all  seasons  indicating 
continuous  redistribution  of  fine  and  coarse 
particles  under  the  influence  of  meteorological 
parameters 

India  1 .  The  mean  and  standard  values  of  the  Angstrom 

exponent,  a  (turbidity  coefficient,  |3),  are  found  to 
be  1 .06  ±  0.33  (0.1 7  ±  0.05)  during  winter, 
0.97±0.35  (0.22  ±  0.04)  during  summer  and 
0.85  ±  0.26  (0.1 6  ±  0.06)  during  monsoon 
2.  The  largest  mean  monthly  AOD  in  February  to 
April  are  associated  with  large  a  values  indicative 
of  fine-mode  particles 
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Table  3 

Studies  of  atmospheric  turbidity  in  America  using  Angstrom  formula. 

Author  Year  Method  Area 

Pinker  et  al.  [83]  2004  Used  an  automatic  sun-tracking  sky  United  States 

radiometer  of  the  CIMEL  type  (Model 
CE  318)  to  measure  both  sun  and  sky 
radiance  at  eight  spectral  channels. 

And  then  used  Dubovik  and  King 
method  (2000)  to  calculate  the  volume 
size  distribution  of  r  from  the  direct 
solar  and  diffuse  sky  radiance 
measurements,  water  vapor 
absorption  band  (940  nm)  to  retrieve 
precipitable  water  vapor  and  the 
Angstrom  exponent  to  calculate 
aerosol  optical  thickness 


Sapkota  and  Dhaubhadel  [84]  2002  Utilized  a  Volz  multispectral  Nepal 

sunphotometer,  equipped  with 
interference  filters  of  3%  bandwidth 
centered  at  440,  500,  640  and  880  nm 
wavelengths  for  the  determination  of 
the  atmospheric  turbidity  over 
Kathmandu  Valley  during  the 
pre-monsoon  period  during  the 
pre-monsoon  period  of  1999  The 
Meteorological  Research  Inc.  (MRI) 
integrating  nephelometer  model  1550 
was  also  operated  to  measure  the  air 
quality  of  the  valley 
The  nephelometer  and  sunphotometer 
were  operated  at  a  height  of  about  45  ft 
from  the  ground  level 

Fox  [58]  1994  Used  the  method  of  Louche  et  al.  [31  ]  Alaska 

to  determine  Angstrom’s  turbidity 
coefficient,  A,  from  measurements  of 
direct  normal  irradiance  (broadband) 
collected  by  the  University  of  Alaska 
Geophysical  Institute  from  1  June  1979 
to  30  June  1983 


Salazar  [85]  2011  Used  the  mean  monthly  global  solar  Argentina 

irradiation  values  measured  at  ten  sites 
in  Argentina  from  the  former 
Argentinean  radiometric  network 
REDSOL  to  estimate  the  monthly 
average  Angstrom  turbidity  coefficient 
|3  values 


Results  and  conclusion 

1.  Annual  means  of  spectral  r  were 
found  to  be  low  over  this  semiarid 
region  and  they  agree  well  with 
historic  data  collected  two  decades  ago 
at  a  nearby  location 

2.  Significant  month-to-month 
variations  in  volume  size  distribution, 
and,  in  particular,  considerable 
variation  in  accumulation  mode 
particles  was  evident 

3.  Monthly  means  of  the  Angstrom 
exponent  were  found  to  be  higher 
during  summer,  suggesting  an 
abundance  of  submicron-sized 
particles 

4.  Precipitable  water  vapor  retrieved 
from  the  940-nm  channel  showed  a 
minimum  of  ~0.4cm  in  February  and  a 
maximum  of  2.5  cm  in  August.  These 
values  were  compared  with 
NOAA/NCEP  Eta  model  forecast  values 
and  found  to  be  in  good  agreement 

1.  About  76.8%  of  the  observed 
Angstrom  coefficient  A  lied  above  0.2 
and  about  14.28%  lied  above  0.4 

2.  About  86.9%  of  the  observed 
scattering  coefficient  bscat  lied  above 
0.21 0  km-1  which  indicates  heavy 
pollution 

3.  The  classes  of  maximum  occurrence 
for  A  and  £>scat  were  found  to  be 
0.22-0.28  and  0.20-0.30  km"1, 
respectively 

4.  The  valley  is  at  the  verge  of 
extremely  heavy  air  pollution 


1.  The  method  was  sensitive  enough  to 
detect  a  seasonal  pattern  along  with 
annual  deviations  associated  with 
hemispheric  scale  disturbance  created 
by  the  El  Chichon  volcanic  eruptions  of 
1982 

2.  Knowledge  of  the  magnitude  and 
variability  of  Angstrom’s  turbidity 
coefficient  will  prove  useful  in 
modeling  solar  radiation  input  to 
northern  ecosystems 

The  monthly  average  value  estimates 
of  turbidity  coefficient  |3  show 
expected  behaviour,  with  values 
varying  according  to  the  season  of  the 
year 


(clean  dry)  atmosphere,  from  the  relative  air  mass.  Abdelrahman 
et  al.  [97]  calculated  the  Linke  factor  and  the  Angstrom  turbidity 
coefficient  f$  at  Dhahran  (Saudi  Arabia).  They  found  that  their  val¬ 
ues  are  larger  in  summer  than  in  winter.  They  constructed  a  linear 
relation  between  the  two  parameters.  Their  coefficients  were  quite 
different  from  those  calculated  at  Avigon  (France)  and  Potsdam 
(West  Germany).  They  attributed  these  differences  to  the  desert 
climate  at  Dhahran  as  stated  by  EI-Hussainy  and  Omran  [39]. 

Grenier  et  al.  [34]  repeated  the  calculation  of  new  <5R(m)  values 
and  proposed  a  fourth  order  polynomial  for  the  function  <5R(m)  sim¬ 
ilar  to  Louche  et  al.  but  with  slightly  different  coefficients,  because 
they  confined  the  range  of  validity  to  m  <  6  in  order  to  achieve  a 
better  fit  in  this  range.  Further,  Grenier  et  al.  separated  the  influ¬ 
ence  of  the  Angstrom  turbidity  coefficient  (3  and  of  the  columnar 
water  vapor  content  w  of  the  atmosphere  on  the  Linke  turbidity 


factor  TL(m)  and  displayed  the  variation  of  TL(m)  with  relative  opti¬ 
cal  air  mass  m  for  various  values  of  /3  and  w  as  parameters.  Grenier 
et  al.  derived  algorithms  for  converting  TL(m)  into  a  standardized 
Linke  turbidity  factor  TL  at  relative  optical  air  mass  m  =  2  with  the 
help  of  polynomials  of  the  third  order  in  TL(m)  whose  coefficients 
are  polynomials  of  the  fourth  order  in  m  as  reported  by  Kasten  [65]. 
Molineaux  et  al.  [90]  reported  that  Grenier  et  al.  [34]  in  the  same 
approach,  added  some  minor  changes  to  the  spectral  absorption 
and  scattering  equations  and  obtained  the  polynomial,  yielding 
very  similar  values  to  Louche’s  [94]  relation.  As  reported  by  Grenier 
et  al.  [95]  recently,  a  model  of  Linke’s  turbidity  factor  was  developed 
by  the  authors  (Grenier  et  al.  [34]),  in  terms  of  spectral  extraterres¬ 
trial  solar  irradiations  and  extinction  coefficients  of  atmospheric 
absorbers.  With  the  help  of  the  model,  the  authors  calculated  the 
theoretical  variations  of  TL  from  the  values  of  the  two  atmospheric 
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Studies  of  atmospheric  turbidity  in  Africa  using  Angstrom  formula. 
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Author 


Year  Method 


Area  Results  and  conclusion 


Tadros  et  al.  [86] 


EI-Hussainy  and  Omran  [39] 


Maduekwe  and  Chendo  [87] 


Trabelsi  and  Masmoudi  [88] 


2002  Estimated  the  ATC,  a  and  A  from  the  r  using 
spectral  broadband  data.  The  data  were 
obtained  from  pyrheliometric  measurements 
in  the  period  1991-96  for  two  sites  of  different 
climatological  and  environmental  view,  Cairo 
in  low  Egypt  and  Aswan  in  Upper  Egypt 
A  technique  used  was  similar  to  that  of  Volz 
(with  two  large  pairs  of  spectral  broadbands 
(bands  1/3  and  bands  2/4)  instead  of  narrow 
bandwidth),  was  applied  to  SMARTS2  and 
SPECTRAL2  models  to  select  the  suitable 
spectral  broadband  for  estimating  ATC 


1998  Calculated  and  analyzed  the  Linke  turbidity  Egypt 

factor,  aerosol  optical  thickness  and  the 
Angstrom  turbidity  parameters  using 
measurements  of  both  total  spectrum  and 
broad  spectral  bands  of  direct  irradiance. 

Relations  were  also  constructed  to  estimate 
both  aerosol  optical  thickness  and  the 
Angstrom  turbidity  coefficient  from  Linke 
turbidity  factor 


1997  Studied  the  trend  in  the  atmospheric  turbidity 

in  Lagos,  a  tropical  city  in  Nigeria,  for  17 
months  between  1990  and  1991  with  a  Volz 
sunphotometer  using  the  Angstrom  turbidity 
coefficient,  /3,  at  500  nm  wavelength 


201 1  Used  the  direct  measurements  taken  by 

pyrheliometer  during  a  year:  July  (2008)  to 
June  (2009)  to  derive  Angstrom  coefficient  (|3) 


Egypt  1.  The  turbidity  A  started  to  increase  from  March, 

and  peak  in  May  with  value  0.31 7  for  Cairo  and 
0.179  for  Aswan.  This  was  attributed  to  Khamsin 
depressions  coming  from  Great  Sahara  associated 
with  advection  of  desert  air  masses  with  an 
extremely  high  concentration  of  aerosols 

2.  The  maximum  of  A  at  summer  months  occurred 
on  June  (0.32)  and  August  (0.1 78)  at  Cairo  and 
Aswan,  respectively.  These  maxima  were  due  to 
hot  and  dry  weather,  wind  dispersion  of  aerosols 
and  dust  raised  from  ground 

3.  The  Sudan  monsoon  trough  prevailing  during 
the  autumn  months  cause  maximum  0.328  at 
September  in  Cairo  and  0.156  at  October  in  Aswan 

4.  The  atmospheric  turbidity  depends  on  the  local 
weather  conditions,  and  the  prevailing  winds  may 
transport  aerosol  particles  from  distant  sources 
depending  on  the  speed  of  wind  and  its  direction. 
The  increase  in  turbidity  during  1992  is  attributed 
to  the  eruption  of  Mount  Pinatubo  in  the 
Philippines  on  15  June  1991 

1.  The  annual  variation  showed  lowest  value  of 
turbidity  in  winter  and  highest  values  in  both 
spring  due  to  the  Khamsin  disturbances  with 
sandstorms,  and  summer  due  to  hot  air  mass  and 
large  water  vapor  contents 

2.  Aerosol  extinction  coefficients  decrease  with  the 
increase  of  both  wavelength  and  optical  air  mass 

3.  Aerosol  optical  thickness  could  be  estimated 
from  the  Linke  turbidity  factor  values  with  errors 
below  20%.  Aerosol  concentration  parameter  A  can 
be  estimated  with  errors  above  50% 

Nigeria  1 .  The  variation  of  /3  indicated  that  the  atmosphere 

was  mostly  turbid  throughout  the  year 

2.  The  aerosols  influencing  the  city  were  both  of 
maritime  and  Saharan  origin  coupled  with  locally 
produced  particulates  from  industries  and 
automobile  combustion 

3.  The  month  of  January  experienced  the  highest 
turbidity,  with  a  mean  value  of  0.497,  while 
November  experienced  the  lowest  aerosol  loading 
on  average  with  the  value  of  0.225 

Tunisia  1 .  The  level  of  turbidity  found  depends  greatly 

upon  the  emissions  and  climate  system  of  each 
area.  The  results  show  the  highest  turbidity  in  the 
summer  season  and  the  low  turbidity  in  the  winter 
season  for  each  index 

2.  Values  of  |3  are  higher  in  the  coastal  region  than 
those  in  the  island  for  all  seasons,  except  in  winter, 
|3  is  almost  the  same  in  both  areas.  The  low  values 
of  |3  in  the  island  of  Kerkennah  can  be  explained  by 
low  amount  of  continental  aerosols  because  in  this 
island  we  have  the  major  party  of  the  territory  is 
occupied  by  lagoons  (Sebkha) 

3.  The  turbidity  of  the  atmosphere  has  a  diurnal 
variation:  it  increases  in  the  afternoon.  This 
situation  can  be  explained  by  a  supplementary 
load  of  the  atmosphere  in  water  vapor  due  to  a  sea 
breeze  circulation  and  the  transport  of  aerosols 
from  local  transport  into  the  site,  essentially 
during  summer  afternoons 


constituents  at  the  origin  of  the  turbidity,  the  water  vapor,  rep¬ 
resented  by  its  content  w,  and  amount  of  aerosol,  represented  by 
Angstrom’s  turbidity  coefficient  /3  and  from  the  relative  optical  air 
mass  m.  Using  the  results  of  the  model,  they  proposed  to  introduce 
a  new  quantity  strictly  representative  of  the  atmospheric  turbidity, 
the  standardized  Linke’s  turbidity  factor,  that  is,  the  usual  factor  TL 
from  the  dependence  on  m  has  been  removed.  Diabate  et  al.  [26] 
reported  that  ESRA  [98]  discusssed  an  approach  for  calculating  TLm, 


and  it  is  found  that  the  analytical  inversion  of  the  clear-sky  model 
is  not  required.  So,  it  is  easier  to  be  implemented  and  operated 
than  that  of  Cucumo  et  al.  [99].  Ineichen  and  Perez  [93]  reported 
Molineaux  et  al.  [90]  noted  that  Louche  and  Grenier’s  expressions 
for  SR(m)  become  divergent,  respectively  for  air  mass  greater 
than  20  and  7.  They  adapted  the  coefficients  of  Linke’s  original 
expression  to  take  into  account  the  absorption  by  the  permanent 
gases.  Cucumo  et  al.  [99]  propose  a  method  for  calculating  TLm  from 
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Table  5 

Number  of  measurements  and  studies  of  atmospheric  turbidity  in  different  areas  in  the  world  using  Linke  turbidity  factor. 

Author  Year  Method  used  Area 

Mavromatakis  and  2007  Used  more  than  a  year’s  data  to  Greece 

Franghiadakis  [102]  determine  the  Linke  turbidity 

coefficient  through  measurements 
of  the  global  and  diffuse  irradiance 
at  5  min  intervals  with  two  CM1 1 
pyranometers  by  Kipp  Zonen  at  an 
altitude  of  122  m 


Chaabane  et  al.  [25]  2004  Used  pyrheliometric  Tunisia 

measurements  in  a  coastal  tourist 
location  in  Tunisia  (Sidi  Bou  Said), 
during  three  summer  months 
(June,  July  and  August  1999)  to 
compute  the  Linke  turbidity  factor 


Chaiwiwatworakul,  2004 

Chirarattananon  [103] 


Diabate  et  al.  [26]  2003 


Derived  the  turbidity  indices, 
namely,  Linke  factor,  Angstrom 
coefficient,  and  illuminance 
turbidity  factor  from  the  record  of 
two  and  half-year  measurement 
derived  directly  from 
measurements  taken  by 
pyrheliometer,  Volz 
sunphotometer  and  beam 
illuminance  meter  at  the  Asian 
Institute  of  Technology  (AIT),  to 
investigate  atmospheric  turbidity 
for  tropical  Thai  sky 


Used  Aguiar  method  1995,  for 
calculating  TL  from  a  time  series  of 
daily  irradiation 


Thailand 


9  stations  in  Egypt,  2  in 
Mozambique  and  Zimbabwe,  1  in 
Algeria,  Tunisia  and  Zambia 


Rapti[104]  2000  Used  global  solar  and  diffuse  sky  Greece 

radiation  data,  measured  with  a 
Moll-Gorczynski  pyranometer  for 
three  values  of  relative 
atmospheric  mass  which 
correspond  to  the  solar  elevations 
of  20°,  28°  and  40°  in  a  coastal 
middle  latitude  Mediterranean 
location,  for  the  time  period 
1973-1976  measured  in  Greece  to 
investigate  the  diurnal  and 
seasonal  variations  of  the 
atmospheric  turbidity  using  Linke’s 
factor 


Results 

1.  Typical  values  of  the  Linke 
turbidity  lied  in  the  range  of 
2.3-3. 5,  while  the  minimum  and 
maximum  measured  values  were 
2.0  and  4.0,  respectively 

2.  Strong  seasonal  effects  on  the 
monthly  averaged  TL  values  were 
not  observed,  while  their  profile 
displayed  a  smooth  transition  from 
lower  to  higher  values 

1.  The  increase  of  TL  is  due  to 
industrial,  increase  of  traffic  in  the 
site  and  transport  of  both  moisture 
from  maritime  air  masses  and 
Saharan  dust  carried  by  sirocco 
winds 

2.  The  atmospheric  turbidity 
depends  on  the  local  weather 
conditions 

3.  Prevailing  winds  may  transport 
aerosol  particles  from  various 
sources  depending  on  the  speed  of 
the  wind  and  its  direction 

1.  Values  of  the  turbidity  of  the 
atmosphere  are  found  to  vary  with 
seasons,  months  and  with  time  of 
day 

2.  Atmospheric  turbidity  is  low  and 
quite  stable  during  dryer  months 
and  increases  in  wet  season  from 
March  to  August 

3.  The  turbidity  increases  from 
morning  to  noontime  and  then 
decreases  in  the  afternoon 

4.  The  annual  mean  values  of 
Angstrom  coefficient,  wavelength 
exponent,  Linke  factor  and 
illuminance  turbidity  factor  are 
0.098, 1.272,  3.306  and  2.818, 
respectively 

1.  TL  is  almost  constant  throughout 
the  year  close  to  the  Mediterranean 
basin  with  values  around  3.5 

2.  Stations  located  in  the 
sub-tropical  Southeastern  part 
exhibit  large  variations  of  TL 

The  range  and  the  behavior  of  the 
monthly  average  values  of  TLm 
throughout  the  year,  demonstrate 
that  TLm  is  related  to  the  climate, 
an  unsurprising  result 

1.  Atmospheric  transparency 
decreased  with  decreasing  relative 
optical  air  mass.  This  leads  to  a 
virtual  increase  of  Linke’s  turbidity 
factor  with  increasing  solar 
elevation 

2.  Diurnal  and  seasonal  variations 
of  atmospheric  transparency  were 
found  with  a  summer  afternoon 
minimum  and  a  winter  morning 
maximum 

3.  Real  diurnal  and  seasonal 
variations  of  atmospheric  turbidity 
were  found  with  a  winter  morning 
minimum  and  a  summer  afternoon 
maximum 

4.  The  largest  diurnal  turbidity 
variation  is  observed  in  October 
(about  6%).  The  annual  variation  of 
turbidity  is  estimated  to  be  about 
40%,  between  the  morning  value  of 
turbidity  in  February  and  the 
afternoon  value  in  June 
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Table  5  ( Continued ) 


Author 

Lopez  and  Batlles  [9] 


Coste  and  Eftimie  [105] 


Kryza  et  al.  [106] 


Year  Method  used 


Area 


Results 


2004  Compared  and  evaluated  the 

performance  of  three  turbidity 
algorithms  -  Dogniaux’s  algorithm, 
Louche’s  algorithm  and 
Gueymard’s  algorithm  -  based  on 
broadband  solar  irradiance 
measurements  in  the  estimation  of 

P 

The  evaluation  of  the  performance 
of  the  models  was  undertaken  by 
graphical  and  statistical  (root  mean 
square  errors  and  mean  bias 
errors)  means 


2010  Use  (Kasten  and  Young  Model,  and 

Remund  and  Page  Model)  to 
calculate  the  Rayleigh  optical 
thickness  (<5r)  and  the  relative 
optical  air  mass  (m)  and  then  to 
determine  the  Linke  turbidity 
factor 

2010  Use  the  formula  proposed  by 

Dogniaux  (1984)  to  determine  the 
Linke  turbidity  factor 


Almeria,  Bondville,  Desert  Rock, 
Fort  Peck,  Goodwin  Creek,  Granada, 
Penn  State,  Table  Mountain 


Romania 


Poland 


1.  Dogniaux’s  and  Louche’s 
algorithms  provide  hourly 
turbidity  estimates  closer  to  each 
other  than  by  using  Gueymard’s 
method 

2.  Differences  on  turbidity 
estimates  by  Dogniaux’s  and 
Louche’s  algorithms  were  notably 
reduced  by  introducing  a 
correction  term  depending  on 
precipitable  water 

3.  Ueymard’s  method  is  the  more 
accurate  and  reliable,  if 
information  about  the  wavelength 
exponent  is  not  available 

4.  Turbidity  algorithms  by 
Dogniaux  [86]  and  Louche  et  al. 

[31  ]  have  shown  to  be  less 
sensitive  to  broadband  solar 
radiation  data  affected  by  cloud 
interference  than  that  by 
Gueymard  and  Vignola  [32] 

1 .  The  turbidity  factor  had  a 
decreasing  variation  during  a  day 

2.  At  the  sunset  time,  the  minimum 
values  of  the  turbidity  factor  were 
obtained 

3.  During  the  morning,  the 
turbidity  factor  has  higher  values 
The  formula  proposed  by  Dogniaux 
(1984)  was  found  to  be  in  good 
agreement  with  measurements  by 
Jacovides  (1997) 


a  time  series  of  daily  irradiation.  This  method  is  based  on  (i)  the 
computation  of  the  direct  and  diffuse  components  by  the  means  of 
the  Collares-Pereira  and  Rabl  correlation,  (ii)  the  modeling  of  each 
component  as  an  analytical  function  of  TLm  and  (iii)  the  inversion  of 
these  models.  The  method  was  applied  to  several  sites  in  Italy.  For 
these  sites,  the  authors  used  the  maximum  daily  irradiation  values 
read  per  month  in  the  ESRA  [100].  These  maxima  are  those  found  in 


the  time-series  of  daily  irradiation  spanning  the  period  1966-1975. 
Additionally,  Walkenhorst  et  al.  [101]  processed  time-series  of 
hourly  means  of  beam  irradiance.  The  Linke  turbidity  factor  was 
calculated  for  the  4  h  around  noon  for  each  day  by  inverting  the 
ESRA  formula  for  the  clear  sky  beam  irradiance.  From  the  resulting 
calculated  set  of  TL  for  each  month,  the  smallest  three  are  averaged 
and  this  mean  serves  as  an  approximation  for  TLm  as  reported  by 


Table  6 

Comparison  between  the  Linke  turbidity  factor  and  Angstrom  turbidity  coefficients. 

Linke  turbidity  factor  TL  The  Angstrom  turbidity  coefficients,  /3  and  a 


1.  Useful  parameter  for  comparison  of  cloudless  atmospheric  conditions 


2.  Refers  to  the  whole  spectrum,  that  is,  overall  spectrally  integrated 
attenuation,  which  includes  presence  of  gaseous  water  vapor  and  aerosols 

3.  It  is  a  key  input  to  several  models  assessing  the  downwelling  irradiance 
under  clear  skies  that  are  used  by  several  communities  in  the  fields  of 
renewable  energies,  climatology,  agro-meteorology,  and  atmospheric 
pollution 

4.  Linke’s  turbidity  factor  has  been  less  used  than  Angstrom  coefficient, 
although  its  experimental  determination  seems  to  be  easier 

5.  Does  not  exclusively  explain  the  turbidity  caused  by  aerosol  because  the 
quantity  is  also  affected  by  the  absorption  of  solar  radiation  in  the  visible 
and  near  IR  regions  by  water  vapor 

6.  Is  not  a  pure  turbidity  coefficient.  It  also  depends  on  air  mass,  m,  and 
precipitable  water,  w  even  when  the  atmospheric  conditions  remain 
constant. 


7.  Intricate  relation  combining  aerosol  and  water  vapor  transmittances 
integrated  over  the  solar  spectrum 

8.  Difficult  to  determine  accurately  and  fails  to  distinguish  the  effects  of  water 
vapor  from  the  absorption  due  to  aerosols 


1.  Angstrom’s  turbidity  coefficient  is  calculated  from  spectral  direct  irradiance 
at  two  wavelengths,  usually  380  and  500  nm,  in  a  part  of  the  solar  spectrum 
where  absorption  is  negligible 

2.  Angstrom’s  turbidity  coefficient  /I  is  usually  determined  under  the 
assumption  of  a  constant  wavelength  exponent  a  =  1 .3 

3.  The  expression  does  not  permit  separate  estimation  of  the  scattering  and 
absorption  optical  depths.  This  causes  important  errors  in  the  estimation  of 
the  diffuse  solar  radiation 

4.  The  determination  of  Angstrom  coefficients  is  limited  to  small  optical  depth 
values.  Therefore  large  optical  depths  that  occur  for  large  zenith  angles  cannot 
be  considered 

5.  The  value  of  /3  can  be  affected  by  a  change  in  the  amount  of  atmospheric 
water  vapor 

6.  The  Angstrom  formula  is  only  a  convenient  approximation.  It  is  not 
necessarily  valid  over  an  extended  spectral  range,  and  Angstrom  parameters 
that  work  properly  in  a  given  range  might  not  be  so  appropriate  in  another. 
Nevertheless,  the  Angstrom  formulas  seem  to  provide  a  good  parameterization 
of  spectral  aerosol  optical  thickness  for  the  400-670  nm  band 

7.  The  Angstrom  method  suffers  from  the  same  problem  that  the  methods 
based  on  photometer  measurements  employing  the  Volz  technique  do,  that  is, 
it  is  highly  dependent  on  the  selection  of  the  spectral  bands  because  different 
band  combinations  produce  different  a  and  /3  values 

8.  The  major  drawback  to  Angstrom’s  formalism  lies  in  the  difficulty  in  getting 
a  precise  evaluation  of  a  and  /I 
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Diabate  et  al.  [26].  Ineichen  and  Perez  [93]  developed  a  new  for¬ 
mulation  for  the  Linke  turbidity  factor  (Tup)  with  the  objective  of 
removing  its  dependence  upon  solar  geometry  and  fully  compati¬ 
ble  with  the  original  formulation  at  air  mass  2.  They  also  developed 
two  new  simple  clear  sky  models  for  global  and  direct  normal 
irradiance. 

Table  5  summarizes  a  number  of  measurements  and  studies  of 
atmospheric  turbidity  in  different  areas  in  the  world  using  Linke 
turbidity  factor. 

4.  Comparison  between  Angstrom’s  turbidity  parameters 
and  Linke  turbidity  coefficient 

Table  6  shows  the  comparison  between  Linke  turbidity  factor  TL 
and  Angstrom  turbidity  coefficients,  /3  and  a. 

5.  Conclusions 

The  results  of  this  overview  indicate: 

1.  The  atmospheric  turbidity  depends  on  local  weather  conditions 
and  on  the  climate  of  the  site.  The  prevailing  winds,  which  may 
transport  moisture  or  aerosol  particles  from  distant  sources,  play 
a  major  role  on  the  seasonal  variation  of  turbidity. 

2.  During  the  last  century,  a  number  of  atmospheric  turbidity 
indices  have  been  introduced  and  several  methods  were  devel¬ 
oped  to  determine  their  values.  The  most  currently  used  are  the 

Angstrom  turbidity  coefficient  /3,  the  Linke  turbidity  factor  IL, 
Schuepp’s  turbidity  coefficient  B  and  Unsworth-Monteith  TU. 

3.  Linke’s  turbidity  factor  TL  represents  turbidity  caused  by  aerosol 
and  by  water  vapor.  TL  is  usually  deduced  from  measurements 
of  the  direct  component  of  solar  radiation  using  pyrheliometric 
methods  during  very  clear  sky  periods,  but  this  kind  of  experi¬ 
mental  data  is  rarely  available,  thus  TL  is  generally  an  estimated 
parameter. 

4.  In  a  pure  Rayleigh  atmosphere  the  turbidity  factor  is  equal  to 
one  (TL  =  1).  The  closest  value  to  this  ideal  value  is  achieved  in 
extremely  clear,  cold  air  at  high  latitudes  (TL  =  2).  However,  for  a 
heavy  polluted  atmosphere  the  turbidity  factor  can  increase  to 
8.  Monthly  mean  values  of  TL  are  normally  about  3  for  temper¬ 
ate  rural  climate,  about  4  between  the  tropics  and  5  for  humid 
polluted  atmospheres.  The  value  for  Linke  turbidity  factor  specif¬ 
ically  recommended  by  the  branch  literature  for  Europe  is  3. 
But,  at  a  local  level  the  value  of  Linke  turbidity  factor  can  vary 
with  the  geographical  coordinates,  meteorological  parameters 
and  climatic  conditions. 

5.  The  Dogniaux  formula  allows  the  calculation  of  the  Linke  tur¬ 
bidity  factor  with  average  errors  of  about  20%  and  quadratic 
errors  of  about  30%  assuming  a  suitable  constant  value  of  the 
Angstrom  coefficient.  As  suggested  by  Grenier  [34,35]  Linke’s 
turbidity  factor  T2. oo  values,  which  correspond  to  the  solar  ele¬ 
vation  of  30°,  can  be  assumed  to  be  independent  of  relative 
atmospheric  mass  and  strictly  representative  of  atmospheric 
turbidity. 

6.  The  main  problem  with  TL  is  that  it  is  difficult  to  be  determined 
accurately  and  fails  to  distinguish  the  effect  of  water  vapor  from 
the  absorption  due  to  aerosol.  Also  it  is  not  truly  independent  of 
the  optical  air  mass  and  varies  diurnally  even  with  an  unchang¬ 
ing  atmosphere.  Moreover,  it  is  subject  to  virtual  daily  variations 
like  broadband  aerosol  optical  depth  <$a  due  to  the  combination 
of  different  parasitic  effects.  Virtual  variations  of  TL  due  to  solar 
zenith  angle,  Z,  alone  are  smoother  than  those  of  <$a  for  Z<  80°  but 
become  steeper  for  larger  zenith  angles.  To  remove  a  part  of  TL’s 
parasitic  variation,  it  is  advisable  to  calculate  it  for  a  reference 
air  mass  value  of  2. 


7.  The  Angstrom  exponent  a  is  very  popular  because  of  the  simplic¬ 
ity  of  the  respective  equation,  because  it  enables  to  interpolate 
or  to  extrapolate  aerosol  optical  properties.  The  Angstrom  tur¬ 
bidity  coefficient  (3  has  long  been  used  as  a  turbidity  index 
and  can  be  determined  by  different  methods  from  spectral 
and  broadband  radiation  measurements  in  the  large  spectrum 
band  not  affected  by  the  water  vapor  absorption  (280<Z> 
630  nm). 

8.  cannot  be  evaluated  unless  a  is  known  a  priori.  This  uncertainty 
in  a  makes  the  determination  of  /3  less  precise  and  discriminat¬ 
ing.  /3  can  be  calculated  in  a  simple  way,  using  the  Katz  formula, 
with  an  average  error  of  about  10%  and  an  average  quadratic 
error  of  30%  if  the  experimental  value  of  the  Linke  turbidity  factor 
is  inserted  into  the  Katz  formula. 

9.  Unlike  the  Linke  turbidity  factor,  Angstrom’s  turbidity  coeffi¬ 
cients  may  be  employed  in  calculations  of  spectral  direct  and 
diffuse  solar  irradiance. 
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